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A curious effect is uncovered by calculating the time evolving probability of reflection of a Gaussian 
wave packet from a rectangular potential barrier while it is perturbed by reducing its height. A time 
interval is found during which this probability of reflection is larger ("superarrivals") than in the 
unperturbed case. This nonclassical effect can be explained by requiring a wave function to act as a 
"field" through which an action, induced by the perturbation of the boundary condition, propagates 
at a speed depending upon the rate of reducing the barrier height. 

PACS number(s): 03.65.Bz 



In recent years a number of interesting investigations 
have been reported on the wave-packet dynamics [Q. 
Here we study it from a hitherto unexplored perspec- 
tive. The reflection/transmission probabihties for the 
scattering of wave packets by various obstacles are usu- 
aUy considered from static or unperturbed potential bar- 
riers. GeneraUy, the time-independent (asymptotic) val- 
ues attained after a complete time evolution are calcu- 
lated. Here an interesting effect is pointed out that oc- 
curs during the time evolution. For this we consider the 
dynamics of a wave packet scattered from a barrier while 
its height is reduced to zero before the asymptotic value 
of reflection probability is reached. 

For an unperturbed barrier, the reflection probability 
for an initially localized wave packet ip {x,t — 0) is calcu- 
lated by considering the wave packet as a superposition 
of plane waves and by writing 



l^ol 



\cl,{p)f\R{p)\'dp 



(1) 



where |i? (p)| is the reflection probability corresponding 
to the plane wave component exp{ipx) and 0(p) is the 
Fourier transform of an initial wave packet ^{x,t ~ 0). 
Since a wave packet evolves in time, |-Ro|^ defined by Eq. 
(1) denotes the time-independent value of reflection prob- 
ability pertaining to a wave packet, which is attained in 
the asymptotic limit {too) of the time evolution. Thus 
|i?oP can be expressed in the following form 



ii?nr = 



\ifi {x,too)\^ dx 



(2) 



where ij} (x,too) is an asymtotic form of the wave packet 
attained by evolving from ^->{x, t = Q) and by being scat- 
tered from a rectangular potential barrier of finite height 
and width. Note that xi lies at the left of the initial pro- 
file of the wave packet such that (x ,t — Q)\^ dx is 



negligible. Equivalence between the expressions (1) and 
(2) in the limit of large t (compared to the time taken 
by the packet to get reflected from the barrier) has been 
checked numerically. At any instant before the constant 
value I-RqI^ is attained, the time evolving reflection prob- 
ability in the region —oo<x<xi\s thus given by 
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(3) 



Now, suppose that during the time evolution of this 
wave packet, the barrier is perturbed by reducing its 
height to zero within a very short interval of time that is 
small compared to the time taken by the reflection prob- 
abihty to attain its asymptotic value |-Ro|^- We compute 
the effects of this "sudden" perturbation on |_R(i)|^. The 
salient features are as follows: (a) A finite time interval 
is found during which \Rit)\ shows an enhancement ( 
"superarrivals") in the perturbed case even though the 
barrier height is reduced. This time interval and the 
amount of enhancement depend on the rate at which the 
barrier height is made zero, (b) We show that the phe- 
nomenon of superarrivals is inherently quantum mechan- 
ical by demonstrating that superarrivals disappear in the 
classical treatment of the problem, (c) The origin of su- 
perarrivals may be understood by considering the wave 
function to act as a "field" through which a disturbance 
from the "kick" provided by perturbing the barrier trav- 
els with a definite speed. 

In order to demonstrate the above features we begin 
by writing the initial wave packet (in the units of ?i = 1 
and m = 1/2) 
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which describes a packet of width ctq centered around 
X = xq with its peak moving with a group velocity 
Vg = 2pQ = ^ towards a rectangular potential barrier. 
The point a;o is chosen such that ip{x,t = 0) has a neg- 
ligible overlap with the barrier. For computing |i?(t)|^ 
given by Eq. (3) the time dependent Schrodinger equa- 
tion is solved by using the numerical methods developed 
by Goldberg, Schey and Schwartz In such a treat- 
ment the parameters are chosen in order to ensure that 
the spreading of a packet is negligible. Here we choose 
xo = —0.3, (To = 0.05/-\/2 and po = SOtt. The barrier is 
centered around Xc — with a width taken to be 0.016. 
For such a width, height of the barrier (V) before per- 
turbation is chosen to be V = 2E, where E is the expec- 
tation value of the energy of the wave packet given by 
Po + j'^o'^ ■ This choice ensures that: (1) The reflection 
probability is close to 1 since we are interested only in 
the reflection probability. (2) At the same time V is not 
too large . This ensures that the reduction of the barrier 
height is not too fast. 

is computed according to Eq. (3) by tak- 
ing various values of x' satisfying the condition x' <= 
Xq — 3(To/a/2- The computed evolution of \R{t)\'^ corre- 
sponds to the building up of reflected particles with time. 
More precisely, it means that a detector located within 
the region — oo < a; < a;' measures by register- 

ing the reflected particles arriving in that region up to 
various instants. First, we compute \R{t)f for the wave 
packet scattered from a static barrier V = 2E. The rel- 
evant curve is shown in Figure 1 which tends towards a 
time-independent value which is the stationary state re- 
flection probability |i?o|^ given by Eq.(l), or equivalently 
by Eq.(2). Next, we proceed to study the consequence 
of reducing the barrier height from V — 2E to V — 0. 
The time evolution of in this case is studied by 

varying the ways in which the barrier height is reduced. 

In the specific cases studied, the potential V goes to 
zero linearly within a switching off time e starting at time 
t = tp chosen to be 8 x 10~* (note that numbers denoting 
the various instants are in terms of time steps; for exam- 
ple, t = 8 X 10""* corresponds to 400 time steps). Here 
e ^to,to being the time required for to attain the 

asymptotic value | -Ro I ^ • The short time span e over which 
the perturbation takes place is given by [tp,tp + e]. tp is 
chosen such that at that instant the overlap of the wave 
packet with the barrier is significant. Figure 1 shows the 
evolution of |-R(i)|^ for various values of e. Varying e 
signifies changing the time span over which the barrier 
height goes to zero which in turn means different rates of 
reduction. Figure 1 reveals that 

\Rpitf = \R,it)\^ t<U (5) 

\Rp{tf > \Rs{t)\' U<t<t, (6) 

\Rp{tf <\Rs{t)\^ t>tc (7) 



where tc is the instant when the two curves cross each 
other, and td is the time from which the curve corre- 
sponding to the perturbed case starts deviating from that 
in the unperturbed case. Here tc > td > tp. 
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FIG. 1. The top curve corresponds to the static case and 
reaches value 1 asymptotically. \R{t)\^ for other curves cor- 
respond to various values of e. The curve with the lowest 
asymptotic value corresponds to the smallest value of e cho- 
sen for this set. As one increases e, superarrivals are slowly 
wiped off. 



As the barrier height is made zero, one does not expect 
at any time an increase in the reflected particle flux com- 
pared to that in the unperturbed case. Nevertheless, the 
inequality (6) shows that there is a finite time interval 
At = tc — td during which the probability of finding re- 
flected particles is more (superarrivals) in the perturbed 
case than when the barrier is left unperturbed (see Fig- 
ure 1). A detector placed in the region x < x' would 
therefore register more counts during this time interval 
At even though the barrier height had been reduced to 
zero prior to that. This effect is essentially quantum me- 
chanical. In order to explicitly show this, we consider 
an initial distribution of particles given by the spreads in 
both position and momentum corresponding to a Gaus- 
sian wave packet. Such particles are now assumed to 
obey the classical equation of motion. We solve numer- 
ically the Liouville equation in this same time-varying 
situation in order to obtain the time-dependent number 
density of particles at the detector. Our results are plot- 
ted in Figure 2 where we show that there are no super- 
arrivals for three different values of e. 

In order to have a quantitative measure of superarrivals 
we define the parameter rj given by 
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where the quantities Ip and Is are defined with respect 
to At during which superarrivals occur 

Ip^ [ \Rpit)fdt (9) 

J At 

Is= [ \Rs{t)\^dt (10) 

J At 
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FIG. 4. Here we show that superarrivals diminish by de- 
creasing the width of the barrier. They completely disappear 
for a small enough barrier width. The three curves for the 
perturbed cases correspond to the widths of 0.016, 0.008, and 
0.004 respectively. 
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FIG. 2. The time- varying reflection probability for the clas- 
sical evolution is plotted for the same values of e as in Figure 
1. The absence of superarrivals in this case demonstrates the 
nonclassical nature of this phenomenon. 



We plot the variation of 77 with respect to e for three 
different detector positions in Figure 3. The effect of 
reducing the barrier width on the magnitude of superar- 
rivals is shown in figure 4. 
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FIG. 3. The magnitude of superarrivals r) diminish with 
an increase in e, the time taken for barrier height reduction. 
This behaviour is seen for three different detector positions 
x'= -0.4, -0.5 and -0.6 respectively. 



The results obtained from Figures 1-4 can be summa- 
rized as follows: (a) There exists a finite time interval 
At during which an increase in the refiection probabil- 
ity (superarrivals) occurs for the perturbed cases com- 
pared to the unperturbed situation, (b) Superarrivals 
are inherently nonclassical. (c) The magnitude of super- 
arrivals T] is appreciable only in cases where the wave 
packet has some significant overlap with the barrier dur- 
ing its switching off. Both rj and At (duration of superar- 
rivals) fall off with increasing e. (d) Superarrivals given 
by rj gradually reduce to zero upon decreasing the barrier 
width, while keeping the initial barrier height B fixed, 
(e) Superarrivals persist^ as we have checked, even if the 
detector is placed at a distance x' greater than Su to the 
left of the initial position of the centre of the wave packet 

Next, we consider the question as to how fast the in- 
fiuence of barrier perturbation travels across the wave 
packet (signal velocity Ve)- Note that even in a classical 
theory the information content of a wave packet does not 
always propagate with the group velocity Vg of a wave 
packet which is usually identified with the velocity of 
the peak of a wave packet ||^. Profiles of the quantum 
wave packet are plotted at various instants in Figure 5. 
An incident packet gets distorted upon hitting the time- 
varying barrier. It splits into two pieces, one of which 
moves towards the right (transmitted particles). The rc- 
fiected packet has a secondary peak shifted towards the 
left. It is thus not possible to uniquely define a group 
velocity Vg for the reflected packet in this case. 



3 



The action due to a local perturbation (reduction of 
barrier height) propagates across the wave packet with a 
signal velocity v^. which affects the time evolving reflec- 
tion probability that can be measured at different points. 
Thus a distant observer who records the growth of reflec- 
tion probability becomes aware of the perturbation of 
the barrier (occuring from an instant tp) at the instant 
td when the time varying reflection probability starts de- 
viating from that in the unperturbed case. Then Ve is 
given by 















^-^w / 









FIG. 5. Snapshots of the wave packet are plotted at four 
different instants of time. The initial narrow Gaussian is heav- 
ily distorted upon striking the barrier. It splits up into two, 
with the reflected part possessing a secondary peak shifted 
towards the detector. 



We compute Ve and Vg =< p > /m for a range of 
parameters and plot v^jvg versus e in Figure 6. Both 
Ai (the duration of superarrivals) and (the signal ve- 
locity) decrease with increasing e (or, decreasing rate of 
perturbation). The magnitude of superarrivals (ry) also 
decreases with increasing e (Figure 3). From such be- 
haviours of rj, At and Ve we infer the following explana- 
tion for the origin of superarrivals. The barrier perturba- 
tion imparts a "kick" on the impinging wave packet which 
spits, and a part of it is reflected with a distortion. A fi- 
nite disturbance proportional to this "kick" or the rate of 
perturbation propagates from the reducing barrier to the 
reflected packet, which results in a proportional magni- 
tude of superarrivals rj. Note that information about the 
barrier perturbation reaches the detector at the instant 
td with a velocity which decreases with the decreasing 
magnitude of impulse imparted to a wave packet. These 
results therefore suggest that information about the bar- 



rier perturbation propagates with a definite speed across 
the wave function which plays the role of a "field" . 
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FIG. 6. The upper curve represents a plot of At (duration 
of superarrivals) versus e. The lower curve is a plot of Ve/vg 
versus e. Here the detector position x' — —0.4. 



To sum up, superarrivals stem essentially from the 
objective reality of a wave function acting as a "field" 
which mediates the propagation of a physical disturbance 
induced by the barrier perturbation. This disturbance 
propagates with a speed that varies according to the rate 
at which the barrier is perturbed. Thus the phenomenon 
of superarrivals has a distinct quantum mechanical sig- 
nificance. Its ramifications call for further studies. In 
particular, different types of perturbations may be stud- 
ied to probe the viability of single particle experiments 
1^,^ for demonstrating this effect. 
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